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Abstract 

We  present  a  model  for  a  flexible  beam  with  time  hysteresis  (Boltzmann  ty]>e  viscoelas¬ 
ticity)  damiiing  and  tip  body.  A  computational  method  for  the  estimation  of  the  damping 
parameters  is  de^•cloped,  and  theoretical  convergence/continuous  dependence  results  are 
given.  An  exami^le  is  presented  in  which  experimental  data  is  used,  demonstrating  the  effi¬ 
cacy  of  the  computational  method  and  the  plausibility  of  the  model  for  predicting  response 
in  dampf^d  structures. 
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I.  Introduction 

In  this  paper  we  present  theoretical  results  for  computational  procedures  develoi)ed  to 
estimate  dami')ing  in  distributed  parameter  models  of  flexible  structures.  These  efforts 
are  part  of  continuing  investigations  of  damping  models  for  composite  material  structures 
that  have  been  the  focus  of  our  attention  for  the  past  several  years.  Our  interest  in 
these  models  grew  out  of  questions  related  to  control  and  stabilization  of  such  structures. 
To  solve  problems  related  to  design,  control  and  staliilization  of  large  flexible  structures, 
one  requires  a)  an  accurate  mathematical  model  of  the  dynamics  of  the  structure,  and 
b)  a  method  or  algorithm  for  estimating  the  parameters  in  the  model.  Most  structures, 
including  the  simplest  of  beams,  exhibit  some  type  of  damping  behavior,  and  it  is  important 
for  the  development  of  accurate  models  to  account  for  the  dam])ing  mechanism.  \''arious 
models  have  been  proposed  for  such  damping  mechanisms  (see  [CP]  and  [R]),  each  grounded 
upon  reasonable  physical  principles  and  yet  leading  to  verj'  different  mathematical  models. 
It  is  important  then,  in  tlie  effort  to  develop  accurate  models  for  flexible  structures,  to  study 
the  various  damping  mechanisms  and  their  implications  for  estimation,  stabilization,  and 
control  problems.  Recently,  Banks,  Inman  et  al  ([BWIC])  reported  on  their  study  of 
Kelvin-Voigt  and  viscous  damping  in  modeling  and  parameter  estimation  for  composite 
beams.  In  related  investigations  Banks  and  Ito  (see  [BI])  developed  a  general  framework 
for  parameter  estimation  (including  numerical  convergence  results)  for  a  large  class  of 
beam  models  with  damping  (including  viscous,  Kelvin-Voigt,  spatial-hysteresis).  A  notable 
exception  that  cannot  be  treated  using  results  of  these  studies  is  viscoelastic  damping  of 
Boltzmann  type  (or  “time  hysteresis”).  Our  focus  in  this  note  will  be  on  this  type  of 
damping.  We  discuss  below  theoretical  and  computational  results  for  the  estimation  of 
parameters  in  a  model  (described  in  section  II)  of  a  beam  with  tip  body  and  Boltzmann 
type  damping.  We  illustrate  use  of  the  resulting  procedures  with  data  from  experiments 
involving  a  composite  material  beam  with  attached  tip  boo_. . 
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II.  Mathematical  Formulation 


The  inatheinatical  model  of  interest  to  us  is  the  following  system  for  the  transverse 
vibrations  of  a  cantilevered  Euler-Bernoulli  beam  with  tip  body  and  Bedtzmann  damping; 

pii„{t.x)  + -^  ~  J  +  s,x)ds^  =  f(i,x)  (2.1) 

0  <  X  <  /,  t  >  0, 


»(bn)=  —(/.())  =  0,  />0, 

ax 


d-i. 


.1)  +  (  J  +  +  IeIhii  -  J  +  s,x)(l.'i I  =  /.•( 


d^u, 


a(0,x)  =  Ut{0,x)  =  0,  0  <  T  <  /. 


Here  u{t,x)  is  the  transverse  displacement  at  time  t  and  position  x  along  the  beam.  Also,  p 
is  the  linear  mass  density,  El  is  the  stiffness,  and  /(t,  x)  represents  external  applied  forces. 
The  boundary  conditions  represent  a  beam  clamped  at  x  =  0  and  free  at  x  =  /,  with 
the  tip  body  dynamics  governed  by  (2.2).  We  assume  that  the  tip  body  has  mass  m  and 
moment  of  inertia  J  about  its  center  of  mass  which  is  assumed  to  be  located  at  a  distance 
c  from  the  tip  of  the  beam  along  the  beam’s  axis  (see  [BR]  for  a  more  detailed  discussion  of 
models  for  beams  with  tip  bodies).  The  functions  k{t)  and  h[t)  denote  externally  applied 
moments  (/.')  and  transversally  applied  forces  (/j)  exerted  on  the  tip  body.  The  delay  length 
r  (thus  the  terminology  ‘"time  hysteresis”)  is  positive,  possibly  infinite.  In  this  paper  we 
assume  that  the  kernel  g{s)  (the  damping  kernel)  has  the  form 


{-sY 


—  r  <  s  <  0, 


n  3  are  positive  constants,  and  p  is  fixed  and  satisfies  0  <  p  <  1. 

Let  q  =  {El  .a,  3.  p.m .  J)  G  Q  C  FI*’  denote  the  parameters  of  interest,  where  the 
admissible  parameter  set  Q  is  a  compact  subset  of  R®.  (The  theory  can  be  modified  to 
include  spatially  varving  El,  but  this  is  not  our  focus  here;  see  [BR],  [BC]}.  It  shouid 
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be  noted  that  we  fix  values  for  r  and  p  rather  than  treat  them  as  parameters.  It  is  also 
interesting  to  treat  ?■  and  p  as  ])arameters,  but  this  requires  the  intrc>duction  of  more 
complicated  parameter-deijendent  state  spaces  (see  [BBC]  for  an  example  of  the  use  of 
such  state-spaces).  We  have  carried  out  n\mierical  studies  in  this  direction  but  will  defer 
any  discussion  on  this  aspect  of  these  models  to  a  sul)sequent  paj>er. 

Here  we  show  (similar  in  spirit  to  the  approaclies  in  [BR],  [BGRW])  how  to  reformulate 
(2.1)  -  (2.2)  as  an  abstract  Cauchy  jiroblem,  and  in  this  framework  we  pose  the  problem 
of  estimating  q  as  a  least-squares  fit  to  data.  Having  done  this,  we  will  in  the  next  section 
discuss  the  construction  of  an  approximation  framework  with  which  the  estimation  problem 
can  be  solved  on  a  computer. 

Tf)  proceed,  let  H  denote  the  Hilbert  space  x  1)  with  inner  product 

where  throughout  we  denote  the  standard  i/-'(0, /)  inner  product  by  (•,•)  .,  and  A  is  the 
2  X  2  jjositive  definite  matrix 

m  me 
me  J  -f  mc^ 

Define  the  Hilbert  space  V  —  {(p,^,^)  £  H  :  (f)  £  /f^(0, /),  (;i'(0)  =  0  =  ci''(0),  p  = 
d(/),  s  =  d'(/)}.  Denoting  elements  of  V  by  $  =  (f)),  *1'  = 

we  take  as  the  V  inner  product  =  (cf>,V')^.  It  is  clear  that  V'  is  continuously 

and  densely  imbedded  in  H,  and  if  we  identify  H  with  its  dual  H*  we  have  the  pivot 
space  framework  V  C  H  =  H*  C  C*.  Next,  define  the  bilinear  form  a  on  V  x  F  by 
^)  = 

It  is  clear  that  the  form  a  is  V'-bounded  and  coercive.  That  is,  there  are  constants 
K  >  0,  C]  >  0,  so  that  for  <5,  G  I  "  we  have 

la(d>.gO|  <  A'I'J’I,.  1^1,,  (2.3) 

and 

<7(5>.  ^)  >  .  (2.4j 
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Now  in  the  usual  fashion  (see  [Kr],  [Sh])  we  can  define  an  operator  A  G  L(V,  V)  by 


for  all  'I'  e  V'. 

Next  define  the  space 


W  =  L2(-r,0;  V;go), 


with  norm  determined  by  the  inner  product 


(Pi.r^),,,  =  I  m(s){T,(s),r2{s})^,ds. 


(2.5) 


Here,  go{s)  =  Qoe^“7(  — s)^,  where  oq,  0o  fixed  values  allowed  for  q  in  the  admissible 
parameter  set  Q.  As  a  notational  convenience,  for  q  =  (El,  a,  p,  m,  J)  G  Q,  w'e  mean  by 
g(q)(s)  the  function  g(q)is)  =  ae^^/(—sY .  We  point  out  that  by  virtue  of  the  definition 
and  the  compactness  of  Q,  the  spaces  L2(— r,  0;  F;  g(q))  are  equivalent  for  all  q  G  Q.  Define 
the  operator  D  ;  domD  C  W  — +  W  on  the  domain 


domD=  {r€^'(-r,0;y;po):  HO)  =  O} 
by  DT  =  dV/ds.  Next  define  the  state  space  Z  =  V  y.  H  y.W  with  norm 


$ 

r 


=  ml  +  ml  +  K- 


Also,  define  the  parameter  dependent  state  operator  A{q)  on  the  domain 

.  ^  \  ^  6  F,  r  G  domD 

domAiq)  -  {  \  ^  ]  GZ:  ^  Jl,giq)(s){r{s)  -  ^)ds)  G  H 


bv 


/ 

(l) 

'I' 


-A(f^  +  i  flg(q)(s)(r(s)  -  ^)ds) 

<l>  +  DT 


(2.G) 
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Finally,  if  we  define  a  state  vector  z{t)  bj’’ 


^-(0=  iUif) 

where  U{t)  =  /),  Vr{t,  /),  u{U  •)),  then  (2.1)  -  (2.2)  can  be  reformulated  as  the  following 

equivalent  abstract  Cauchy  problem  on  Z: 

i(t)  =  +  F{i),  z{0)  =  Z,  (2.7) 

where  F{t)  =  ^F(/)^  with  F{t)  =  j  ^  ^ 

poscdness  result  for  (2.7),  we  make  the  following  assumption  on  the  admissible  parameters. 
Al)  There  exists  e  >  0  such  that 


p  >  e,  El 


-/_>(.)( 


s)ds  >  e  for  all  q  £  Q. 


This  is  a  physically  reasonable  assumption  (see  [Wa]),  basically  guaranteeing  that  the 
structure  is  a  solid.  We  shall  have  occasion  to  consider  the  space  Zg,  which  is  the  Hilbert 
space  consisting  of  the  elements  of  Z  equipped  with  the  inner  product 

(($i,^/i,ri),($2,^'2,r2)),  =iEI-  9{q){s)ds)a{^u^2)  +  p{^u^2)„ 

+  j  9{Q){s)(T{Ti{s),T2{s))ds, 

for  q  =  {EI,a,  P,  p,m,  J)  £  Q.  Since  Q  is  compact,  one  can  readily  see  that  these  inner 
products  are  equivalent  on  Z  for  all  q  £  Q.  Next  we  collect  some  results  concerning  the 
operators  A{q)  and  D  which  we  have  defined. 


Theorem  2.1. 

i)  Under  assumption  Al),  the  operator  A{q)  is  dissipative  in  Zg  for  each  q  £  Q.  Further, 
A{q)  generates  a  Co-semigroup  T{q){t)  on  Z  which  satisfies  \T{q)(t)\^  <  Me‘^\  where  the 
constants  M,ui  are  independent  of  q  £  Q. 
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ii)  The  operator  D  is  dissipative  on  Tl’,  and  the  resolvent  (A/  —  D)  ^  exists  for  A  >  0. 


These  results  follow  by  extending  the  arguments  found  in  [FI].  We  can  thus  conclude  that 
the  system  (2.7)  is  well- posed. 

We  may  now  formulate  the  estimation  problem  of  interest  as  follows:  Find  q  £  Q  that 
minimizes 

=  (2.8) 

The  values  /z,  are  observations  of  the  tip  velocity  of  the  beam  (obtained,  say,  via  a  laser 
vibrometer)  at  time  t,,  and  v(q){ti,l)  is  the  fourth  component  of  the  state  vector  z(q)(t) 
(i.e.,  the  first  component  of  where  z{q){t)  is  the  mild  solution  of  (2.7)  given  by 

=  T{q){t)zo  +  f  T{q){t  -  s)F{s)ds.  (2.9) 

Jo 

That  this  inverse  problem  is  well-posed  will  follow  from  the  compactness  of  Q  if  we  show 
that  J  is  continuous  in  q.  But  in  view  of  (2.9),  under  reasonable  assumptions  on  the  initial 
data  Zq  and  the  forcing  term  F{t),  this  follows  if  we  show  that  T{q){t)2  — >  T{q){t)z  for 
arbitrary  q  q  and  z  £  Z .  The  following  theorem  provides  the  desired  results. 

Theorem  2.2.  If  q  q,  then  T{q){t)z  — *  T{q){t)z  for  all  z  £  Z,  uniformly  in  bounded 
t -intervals. 

Proof  Since  the  generators  .^(<7)  are  dissipative  on  Zq,  the  result  follows  from  the  well 
known  Trotter-Kato  theorem  (see  [Pa])  once  we  show  continuity  in  q  of  the  resolvents 
Rx  {A{q)).  To  this  end,  let 

.-  =  .-(j)=  j  =BU.4(g))i  (2.10) 

and 

s  =  Hq)=  =Rx{A{q))x  (2.11) 


We  must  show  that,  for  some  A  >  u;,  z  z  as  q  q.  Fix  A  >  max{(j,',0}.  From  (2.10) 


wc  ha\-e 


A*^  -  ^ 


A^/  +  yA  (^EIA>  +  j  g{q){s){T{s)  -  ^)ds^  =  r; 


AF  -  -  DF  =  C- 


Substituting  the  first  and  third  equations  into  the  second,  we  have 

A"$  +  J  (^EI  -  J  g{q)is)e^\L^  AA> 


=  ^  ^  giq)(s)A  [(A7  -  (C  -  O]  ds. 


(2.12) 


Similarly,  from  (2.11)  we  have 

+  J  (^EI  -  j  g(q){s)e^^ds^  A^ 


=  g  +  \i--.j'  g(q){s)A  [(A/  -  D)-\C  -  O]  ds. 


(2.13) 


Define  the  bilinear  form  <j>  on  F  x  V  by 


aA(?)($,  ^')  =  A-($,  f  9{q){s)t^'ds^  a($,  <l>).  (2.14) 

Since  A  >  0,  it  is  clear  that  ax  is  F-coercive  [this  follows  from  Al)  since  —  g{q){s)ds  < 
—  g{q)is)e^^ ds].  That  is,  there  exists  c  >  0  and  A  so  that 


c)4>)^^  4-  AK>1^^  <  cta(9)(4>,  4>)  for  all  g  6  Q- 
Next  note  that  it  follows  from  the  form  of  g{q){s)  that  for  fi  E  i7^(— r,0), 

J  \9{q)i^)  -  9{q){s)\nis)ds  — ^  0 


(2.15) 


(2.16) 


as  q  —*  q.  From  this  it  follows  that 


\axiq){^.  ^ - 

I  p  p 

rO 


J-r'P  P  '  '  ' 

<(U(q^q)-\^'\,.-\'i'\,,  (2.17) 


where  dx{q.q)  — -  0  as  q  —  q. 

Xow.  since  <!'  —  <!'  =  A(<J>  —  4*),  and  F  —  f  =  (1  —  finished  with  the 

proof  if  we  show  that  |d>  —  — +  0  as  ^  — +  To  do  this,  we  make  use  of  the  ahove-defined 

coercive  form  ax.  From  (2.15)  and  (2.12)  -  (2.13),  we  have 


c|^>  _  4>|^  +  A|d'  -  5'|*^  <  axiq){^'  -  -  ^>) 

=  (7a((7)(<I',^  -  I')  -  <ta(9)($,^  -  ^) 

axiq){^,  $-!>)-  0\{q){^,  ^>  -  $) 


+ 


<  dx{q,q)\^\,A^  - 

-q{q)i^)  -  ^9{q)is) 


+  K 


L 


i: 


,X(a-e) 


[C(^)-^]cZ0  ds\^-^\ 


Cancelling  a  term  1^>  —  ^>1^^  from  each  side  of  the  inequality,  we  obtain  an  estimate  from 
which  the  result  follows. 


Let  us  make  a  few  remarks  about  the  proof.  Many  proofs  of  convergence  of  a  sequence 
of  semigroups  involve  the  application  of  some  version  of  the  Trotter-Kato  theorem.  This 
typically  involves  showing  convergence  either  of  the  generators  of  the  semigroups  or  of  their 
resolvents,  depending  upon  the  version  of  the  Trotter-Kato  theorem  which  is  used.  For  pa¬ 
rameter  estimation  problems  such  as  we  are  considering  here,  the  “resolvent  convergence” 
form  of  the  Trotter-Kato  theorem  is  preferable  as  it  leads  to  convergence  arguments  similar 
in  spirit  to  those  of  the  “variational  approach”  of  finite  element  theory  (i.e.,  the  definition 
and  use  of  the  form  ax)-  Also,  for  problems  involving  spatially  varying  parameters,  this 
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\'c‘rsi()n  allows  one  to  impose  weaker  smoothness  requirements  on  the  parameters  in  order 
to  ju'ove  convergence  results  (for  a  more  detailed  discussion,  see  [B],  [BI]). 

These  remtirks  will  he  relevant  later  in  the  j^aper  when,  after  constructing  an  ai>proxi- 
mtitiou  scheme  for  our  estimtition  ]iroi>lem.  we  again  a])ply  the  Trotter-Kato  theorem  to 
aigue  convergencf'  of  a  s('f[uence  of  a])])roxiniatiug  seniigrou])S. 

III.  AlM'IiOXl.MATION  FUAMEWOUK 

Constructing  a  finite  dimensional  a])i)roxiniation  scheme  for  the  ])ar;uneter  estimtition 
ju'ohh  in  (2.S)  -  (2.9)  iin'olves  aj^iroximating  the  state  operator  A{q)  and  the  state  sjiace 
Z.  (Were  Q  also  infinite  dimensiontil.  which  is  often  the  case  for  estimation  jirohlems, 
we  would  idso  need  finite  dimensionid  sj^aces  to  tq^ij^roximate  Q.  Treatment  of  this 
case  is  iiossihle  hut  we  will  not  pursue  it  here.)  We  observe  that  approximating  A(q)  and 
Z  im'olves  two  stages  -  apjiro.ximation  in  the  spatial  variable,  and  ajipro.ximation  in  the 
dehiy  variable.  Wc  will  see  that  the  two  ajiproximation  sttiges  arc  independent,  so  that 
our  task  in  constructing  a  scheme  for  (2.S)  -  (2.9)  reduces  to  choosing  a  reasonable  spatial 
approximation  (such  as  finite  elements)  and  a  reasonable  delay  approximation  scheme 
(there  are  several  in  the  literature  (see  [BB].  [BK],  [IK])).  That  is,  there  is  no  conditional 
relationship  between  the  spatial  and  delay  variable  approximations.  For  definiteness,  we 
use  the  index  .V  for  the  spatial  variable,  and  M  for  the  delay  variable. 

Let  us  consider  first  the  approximation  in  the  sjiatial  variable.  This  amounts  to  choosing 
a  sequence  of  finite  dimensional  spaces  satisfying  C  ^  '  C  H  which  approximate  the 
spaces  \  '  and  H  (see  .42)  below).  Once  the  are  given,  the  rest  of  the  approximation  in 
the  spatial  variable  follows  naturally.  That  is.  the  form  a  defines  operators  .4’'  :  H’'  — ♦  H’' 
by 

(.4-'  <J>.  =  <7(4>.  ^0  for  all  «>,  <l’  €  H^' . 

.41so.  define  spaces  IF'  C  B’  by  ir'  =  L2i—r.O:H''\go)  «ind  Z''  C  Z  by 
Z''  =  X  H‘^  X  irC  Let  P,)  ;  H  :  V'  ^  and  :  IF  IF’'' 

denote  the  respective  orthogonal  projections,  and  define  the  operator  D’'  on  the  domain 


domP'  =  (domP)  H  11 


bv 


D^'T  =  DT. 

Finally,  define  the  state  operator  >^'''(9)  on  the  domain 

-i>\ 

G  :  r  G  doml?'^’ 

r/ 

1)V 

/4>\  /  <1^ 

-  i  /!,5(9)(.<^)-4"(r(.)  -  i>)ds  )  .  (3.1) 

Vry  \  <i>  +  D'''r  J 

W'e  note  that  each  A^'iq)  generates  a  Co-semigroup  on  Z''  (this  involves  a  slight  modifi¬ 
cation  of  the  proof  that  A{q)  generates  a  Co-semigroup). 

Next,  we  consider  approximation  in  the  delay  variable.  That  is,  we  apply  one  of  several 
schemes  in  the  literature  for  delay  ecjuations  to  the  operators  A^  {q)  and  spaces  Z''' .  This 
amounts  to  defining  (  for  each  A')  finite  dimensional  subspaces  C  IF'^'  and  operators 

D^-^'  :  However,  we  are  interested  in  imposing  conditions  on  the  delay 

variable  approximation  which  are  independent  of  the  spatial  approximation  (hence  inde¬ 
pendent  of  A').  In  order  to  do  this,  we  will  consider  the  generic  space  IF  =  L2(  — r,  0;  R;  go), 
where  R  =  R' ,  and  the  generic  operator  D  defined  on  the  domain 

domA>  =  {w  G  (— r, 0;  R;  </{))  :  u:(0)  =  O] 


dom>l^’(5)  = 


Iw 


-  dw 
Dw  = 

ds 


A  delay  approximation  scheme  consists  of  finite  dimensional  subspaces  iF^'  C  IF  and 
operators  ;  iF*^'  — +  iF'''  satisfying  certain  conditions  (see  A3)  and  A4)  below)  which 
are  independent  of  A'.  However,  although  iF'*'  and  O’"  are  independent  of  A',  for  each  N 
they  define  in  a  natural  way  subspaccs  IF'^’^'  C  IF^'  and  operators  D'""  ''  ;  IF^'^'  — +  IF'^'  ‘'F 
That  is.  suppose 

=  eRj. 
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where  c*'  are  given  basis  functions  for  M’*'.  Then  define  by 

^  (3.2) 

Iff/'  =  |E.«.^r  eRj  ,  where  the  are  basis  functions  for  ,  then  any  element 
r  of  ir'  '''  can  be  written  as 


r  =  E,  r^xj'  =  E,  (E,  «;  €  R. 

Hence  we  define  D’'  by 

”(^  Ej(E.‘-;sr)c;']  =  Ei  {r>‘’{E,oic‘;))B;-.  (3,3) 

In  a  similar  manner,  the  orthogonal  projection  fl*‘  :  W  — >  W'"’  defines  naturally  the 
orthogonal  projection  H'^  '''  :  IT''’  — »  This  completes  the  cipproximation  of  the 

delay  variable, 

Ha\'ing  the  appi'oximation  of  the  spatial  variable  (the  spaces  H''  )  and  the  delay  variable 
(the  spaces  H  '''  and  operators  D^'),  we  define  the  finite  dimensional  state  space  Z'"  ''’  by 

and  state  operator  on  Z'’'  ’'’  by 

(  4'  I  =  (  -“.4''$-i/!,5(5)(.<).V(r(.«)-4'K'  )  .  (3.5) 

\T  J  \  +  / 

Finally,  we  state  the  finite  dimensional  estimation  problem  of  interest:  Find  q  E  Q  that 
minimizes 

E  P  (3.6) 

I 

where  iq){'(-J)  is  the  fourth  component  of  the  state  vector  {q){t)  given  by 

c^-«(ry)(t)  =  T^'^'{q){i)P^''^^Zo  +  /  T^''^’{q)it  -  s)P^'^' F{s)ds.  (3.7) 

Jo 
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Here  T^-^’{q)  is  the  semigroup  generated  })y  (q)  and  ^  is  the  canonical  orthogonal 

projection  of  Z  onto  Z  ’'  Using  arguments  similar  to  those  in  the  proof  of  Theorem  2.2, 
it  is  readily  seen  that  Z'^  and  JZ  "'  are  continuous  in  q  for  q  in  the  compact  set  Q.  Thus, 
for  each  .V,  M .  the  problem  for  (3.G)  -  (3.7)  has  a  solution  q^  .  To  obtain  convergence  and 
continuous  dependence  (with  res])ect  to  the  observations  n, )  of  the  ])arameter  estimates 
obtained  in  the  finite  dimensional  problems  for  (3.G),  it  suffices  under  the  assumption 
that  Q  is  compact  (for  details,  .see  [B],  [BKu])  to  argue:  For  arbitraiy  {9^  '''}  C  Q  with 
q^'^'  — ‘  17,  we  have  z'' {q’' ){i)  — »  z{q){t)  where  and  r  are  given  by  (3.7)  and  (2.9), 
respectively.  It  is  readily  argued  that  this  convergence  holds  under  reasonable  assumiitions 
on  Co  inid  F{t)  if  one  first  argues  that  T''"*  ){1)P2^  z  — *  T{q){i)z  for  arbitrary 

fj''"''  q  in  Q  and  z  E  Z.  We  give  the  arguments  for  this  latter  convergence  in  the  next 
section. 

IV.  Convergence  Analysis 

To  give  the  desired  convergence  arguments,  we  need  several  approximation  assumptions, 
which  we  now  state. 

A2)  C  V  C  H  satisfies  \P^^'  —  4'|  —♦0  and  —  $1  — >  0  as  A"  00. 

_  ''  '' 

A3)  C  W  satisfies  —  tel  — ►  0  as  M  — >  00. 

«■ 

A4)  i)^'  ;  satisfies 

i)  D'’  is  dissipative  on  (uniformly  in  M), 

ii)  ( XI  —  tr  — ►  {XI  —  D)~^w  as  M  00  for  Re  A  >  0. 


Remark  4.1  Recall  that  to  argue  U- coerciveness  of  the  form  a>,  we  used  the  fact  that 
g(^)c^^ds  <  g{s)ds.  The  function  may  be  characterized  by 

c^^  =  f(s)~X{XI-D)-^f{s), 

where  f{s)  =  1.  The  delay  approximation  scheme  defines  a  natural  approximation  of 
given  by 

e^'(s)  =  n''/-  A(A/-  £)^')-'n"7'. 
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In  order  to  argue  I '-coerciveness  of  a  form  which  we  define  below,  we  need  to 

satisfy  an  inequality  analogous  to  that  given  above: 


for  all  sufficiently  large  AI . 


(4.1) 


Before  proving  the  main  semigroup  convergence  result  of  this  section,  we  state  the  fol¬ 
lowing  lemma. 

Lemma  4.1.  Assume  that  A3)  and  A4)  hold,  and  that  and  arc  given  l)y 

(3.2)  and  (3.3).  Then  D’'  "'  is  dissipative  on  (uniformly  in  N.M)  and 


-^o  (4.2) 


as  M  — ♦  cx)  for  each  N  and  each  F  6 


The  proof  is  straightforward.  We  now'  prove  the  semigroup  convergence  result  referred  to 
above. 


Theorem  4.2.  Assume  that  Al)  -  A4)  and  (4.1)  hold.  If  —*  q  in  Q,  then 

jA.Af  (^A.Ar  )pA’.M^  — >  T{q){t)z  for  all  z  E  Z,  uniformly  in  bounded  t-intervals. 


Proof  First,  the  dissipativeness  of  implies  that  the  A^’^*(q)  are  dissipative  on 

Zg,  uniformly  in  N  and  M.  Thus,  w'e  can  apply  the  “resolvent  convergence”  form  of  the 
Trotter-Kato  theorem,  as  in  the  proof  of  Theorem  2.2.  Let 

j  =  R^(A(q))x  (4.3) 


/  A 

.A-.Ar  ^  (4.4) 

where 

/n  /  \  /  r 

X  =  U  and 

\cj  \pr'c/ 
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Solving  (4.3)  and  (4.4)  for  we  find 

A-4>  (/{cj){s)c^-'ds^ 

=  v  +  Xi--^  /° 

and 

-  f  j(?"'")(.')ca(.'M^)  .4"*^ 

=  /”  Sll'’'  *'  )<■')■■»'''  If^-'  -  ®"'’'  -  f")]  (4  ''-) 

As  in  the  proof  of  Theorem  2.2,  it  is  sufficient  to  argue  that  |4>'^  ~  — +  0  as  —*  q. 

To  tliis  end,  define  a  bilinear  form  ct”  on  V  X  V'  by 

=  A"($,'S)„  +UeI-  g(g}{s}c'‘^{s)ds'j  al*,!')-  (4.7) 

From  the  technical  condition  (4.1)  it  follows  that  <7^  is  V-coercive,  uniformly  in  M  and  q, 
for  some  A  >  0.  That  is,  there  exists  C  >  0  so  that 

C|4>|^,  +  A|4>r;  <ar (?)($,$)  (4.8) 

for  all  q  £  Q.  Also,  it  is  easy  to  see  that  for  all  €  V , 

kA(?)('»,  S 4-)!  <  ^  9(q)(s)  [e^'  -  £t'(^)]  ds  •  A-  •  |4.|„  lil,, 

<  A(M)|4'|,,|»i,, 

where  K{M)  — >  0  as  M  — >  oo.  Now,  let  satisfy  — +  0  as  .V  — >  oc 

(existence  of  a  $'''  follows  from  A2)).  Thus  it  is  sufficient  to  argue  that  —  I*''  — +  0. 
Letting  6'^  =  Ave  have 

C|$^'  _  4,n|2  +  ^|$A'  _  |,a'|2  < 

=  <(9--)($-,6")-cr,(g)($,6-)  +  <rA(g)(«J',<5-)-<(g"’*')(e",<^") 
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i 


-f 

P  J-r 


+  -  /  ((A/  -  Dr\(:  -  ds  +  a,{q){^,6n  -  (9"  "  )(^^«^") 


K 


+ 


7-r^ 


g[cr-'^'){s)\ixi-D-^-^‘r\C'^‘-c')-i\i-Dr'ic-i)\,M^^\, 

+  (7A(ry)(3',<5")-aA(9'''^')(4>,,S^-) 

+  )(*J’  - 


(a) 

(b) 

(c) 
(cl) 

(e) 

(f) 


The  proof  is  comijlcted  by  tipplying  A2)  to  («),  (4.2)  to  (b),  (2.1G)  to  (c).  (2.17)  to  (d), 
(4.9)  to  (e),  and  the  uniform  V'-boundedness  of  to  (/). 

V.  Computational  Results 

In  this  section  we  describe  some  of  our  computational  results  for  the  Boltzmann  damping 
model  under  consideration.  Recall  that  for  the  approximation  scheme  which  we  have 
described,  the  underlying  spatial  and  delay  discretizations  should  satisfy  assumptions  A2) 
-  A4)  and  (4.1).  For  most  of  the  numerical  experiments  which  we  have  run,  we  have  based 
our  approximations  on  cubic  splines  for  the  spatial  variable  and  a  version  of  the  so-called 
averaging  scheme  for  the  delay  variable. 

More  specifically,  for  a  beam  of  length  1  and  a  given  positive  integer  A’,  consider  the 
partition  of  the  spatial  interval  [0,/]  given  by  A^'  =  {x^j-Co,  where  x,  =  j:  ■  Let 
S^(A''  )  denote  the  set  of  cubic  splines  with  knots  A^,  and  the  standard  cubic 

B-spline  basis  set  for  5^(A'' )  (see  [Pr]).  Next,  we  can  in  the  usual  way  define  a  new 
set  of  functions  so  that  the  fixed  end  (x  =  0)  boundary  conditions  are  satisfied. 

That  is,  define  —  b^  -1-  26^^,  and  b^  —  6f  for  t  =  2, 3, . . .  ,  A’’  4-  1.  Then  define 

=  (bj'(l),  {l),bj'),  and  set  H''  =  span  {Bf,...  For  this  choice  of 

assumption  A2)  is  satisfied. 
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For  the  uelay  variable,  we  use  the  “non-uniform  mesh”  version  of  the  averaging  scheme 
(see  [FI]).  That  is.  for  a  delay  length  r,  positive  integer  M  and  kernel  function  g{s),  set 
C  =  Then  the  partition 


7  ^  ^  M  -  l  ^ 


<  <  /"  =  0 


is  defined  uniquely  by  the  condition  g(s)ds  =  For  the  kernel  functions  g{s)  under 

} 

consideration,  which  are  increasing  and  have  a  singularity  at  s  =  0,  this  results  in  a  finer 
mesh  near  .c  =  0.  Next,  let  denote  the  characteristic  function  on  the  interval 
Define  the  spare  H’*'  =  span  {x*N...  Finalh',  the  operator  :  11''''  — >  IF'''  is 

defined  bv 


j=i  j=i  D-i  } 


^  j  - 1  A/ 

A»  Aj  ! 


where  Go  =  0.  It  can  be  shown  (see  [FI])  that  IF*'  and  D'''  satisfy  A.3)-A4)  and  that  (4.1) 
is  satisfied. 

We  ha\-e  conducted  a  number  of  numerical  tests  using  simulated  data  in  which  we  ob¬ 
served  the  convergence  and  efficiency  of  the  scheme  for  identifying  the  parameters  El,  q, 
3,  f).  m  and  J.  However,  here  we  present  an  application  using  experimental  data  in  order 
to  demonstrate  that  the  Boltzmann  damping  model  can  provide  a  very  good  prediction  of 
vibration  response  in  certain  flexible  structures.  In  particular,  we  use  data  from  a  series  of 
experiments  carried  out  at  the  Mechanical  Systems  Laboratory  at  SUNV  at  Buffalo.  (We 
gratefully  acknowledge  Dr.  D.J.  Inman,  H.  Cudney,  and  Z.  Liang  for  their  collaboration 
on  these  experiments.)  In  these  experiments,  a  cantilevered  composite  beam  with  tip  body 
was  excited  from  an  initial  rest  configuration  via  an  impulse  from  a  force  hammer.  A  laser 
vibrometer  was  used  to  collect  measui'ements  of  the  tip  velocity  of  the  beam  in  response 
to  the  force.  These  measurements  were  used  as  observations  in  our  least-squares  criterion 
J  (recall  (2.8)). 

For  the  particular  experiment  presented  here,  we  used  a  composite  beam  of  length 
/  =  1.9CS  ft.,  width  0.3345  ft.,  and  thickness  0.01968  ft.,  with  a  tip  body  (a  cylinder 
of  radius  0.08036  ft.  and  height  0.27224  ft.)  with  mass  m  —  0.075714  slug  and  moment  of 
inertia  J  given  by  J  =  cjm  with  Cj  =  0.088265  ft.  The  beam  density  was  p  =  0.021441 
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slug/ft.  We  observed  the  tip  velocity  of  the  beam  after  it  was  hit  with  an  impulse  hammer 
at  Xj,  =  0.9S4  ft.  The  shape,  height  and  duration  of  the  input  were  modelled  using 
the  input  signal  from  a  transducer  in  the  hammer.  We  recorded  2048  time  observations 
at  At  =  0.001953  seconds  intervals,  thus  giving  time  domain  data  over  the  ])eriod  [0.4] 
seconds. 

The  actual  structure  (a  composite  material  beam  with  a  cylindriciil  body  attached  near 
the  free  end)  used  in  our  experiments  entailed  the  “tip”  body  being  attached  with  its 
center  of  mass  a*  1.9352  ft.  from  the  clamped  end  of  the  beam  of  1.9CS  ft.  in  length.  So 
in  fact  we  did  not  have  a  “tip  l)ody"  as  descril>ed  in  the  model  above.  To  compensate  for 
this  inaccuracy  in  the  model,  we  took  c  =  0  but  allowed  Cj  and  m  (i.e.,  J  and  ni)  along 
with  El,  a,  and  3  to  be  estimated  Iw  our  algorithm.  We  would  expect  to  (and  did)  obtain 
“effective”  values  of  Cj  and  in  that  are  slight  perturbations  of  the  measured  values  given 
above. 

We  then  used  data  from  the  experiment  (1018  obser\’ations  in  the  time  period  [0.01172,2] 
-  we  did  not  use  any  observations  from  the  time  period  for  the  hammer  hit)  with  the 
procedure  we  have  described  in  order  to  estimate  the  stiffness  parameter  El,  the  damping 
parameters  a  and  ,3,  and  Cj  and  m.  In  these  computations  we  chose  p  =  0.5  and  r  = 
0.01.  We  obtained  the  values  El  =  103.892  lb  •  ft^,  o  =  48.409,  3  =  654.109.  = 

0.091026  ft.,  and  m  =  0.070044  slug.  A  value  of  El  was  also  determined,  using  standard 
modal  teclmiques,  to  be  107.103  lb  •  ft^.  Thus  our  value  is  well  within  the  10%  range 
which  is  generally  thought  to  be  the  accuracy  of  the  modal  techniques.  The  estimated 
values  for  Cj  and  m  are  also  within  10%  of  the  measured  values.  We  point  out  that  an 
heuristic  physical  interpretation  of  the  large  value  of  /?  is  that  the  material  has  a  “memory” 
which  fades  very-  rapidly.  This  type  of  short  “memory”  makes  sense  for  materials  which 
are  “more  elastic  than  viscous”  as  our  experimental  beam  appeared  to  be.  In  fact,  for 
the  case  of  no  “memory”  at  all  (i.e.  g{s)  =  0),  our  model  reduces  to  that  of  a  purely 
elastic  Euler-Bemoulli  beam.  Thus  the  presence  of  a  small  delay  (the  memory)  in  the 
mathematical  model  drastically  changes  the  qualitative  behavior  of  the  model  (from  elastic 
to  viscoelastic).  (That  small  delays  can  affect  behavior  of  a  model  has  been  observed 
elsewhere;  see  [DLP]). 
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We  depict  in  Figures  1,  la,  and  2  the  time  domain  response  and  FFT  for  our  beam 
model  and  exi)eriment.  In  these  figures,  the  solid  line  graphs  represent  the  experimental 
data  while  the  dashed  lines  represent  the  model  response  corresponding  to  our  estimated 
parameter  values.  It  is  clear  that  we  have  obtained  a  verj’  good  fit  to  data  with  this  model. 
Other  similar  experiments  and  computational  studies  ([BWIC],  [DFWIC])  illustrate  the 
usefulness  of  the  techniques  described  in  this  paper  in  studying  damping  in  distributed 
parameter  models  for  structures. 
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